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1 .  jn  t  r  (^1  'ct  ion 

In  studies  of  wave  propagation  in  magnetized  plasmas  it  is  often  of 
interest  to  excite  low  frequency  modes  from  a  remote  location  under  con¬ 
trollable  conditions.  These  modes  can  be  used  for  a  variety  of  purposes, 

e.g.,  as  a  tool  for  diagnosing  plasma  properties, to  provide  a  seed  for 

(2  3) 

growtli  of  low  frequency  turbulence,^  *  or  perhaps  to  establish  a  low  fre- 
ipiency  communications  link,  as  may  be  the  case  in  the  ionosphere. Regard¬ 
less  of  the  particular  application,  a  common  problem  is  that  at  low  frequencies 
it  is  difficult  to  excite  an  external  electromagnetic  signal  that  penetrates 
deeply  into  the  plasma  where  the  propagation  of  the  low  frequency  mode  is 
desired,  lor  this  reason  one  is  lead  to  consider  processes  in  which  the 
low  frequency  mode  is  indirectly  excited  by  energy  sources  which  penetrate 
the  plasma  more  readily.  Two  obvious  choices  are  particle  streams  and  high 
frequency  waves.  As  is  well  known,  both  of  these  sources  can  trigger  a 
variety  of  parametric  instabilities  whose  decay  channels  can  coincide  with 
the  desired  low  frequency  mode.  An  intrinsic  drawback  of  the  parametric 
excitation  scheme  is  that  the  amplitude  and  phase  of  the  low  frequency  mode 
cannot  be  controlled  by  the  experimentalist,  and  are  difficult  to  predict 
theoretically.  Essentially,  because  the  saturation  of  the  instability  is 
the  result  of  complicated  nonlinear  interactions. 

A  more  controllable  process  consists  of  beating  two  high  frequency 
pump  waves  of  frequency  uij  and  to  drive  a  low  frequency  mode  at  m  - 

resonantly.  In  this  process  the  power  transferred  to  the  low  frequency 
mode  depends  on  the  product  of  the  individual  pump  powers,  and  as  is  illus¬ 
trated  in  this  study,  the  transfer  coefficient  can  be  calculated  from  first 
principles.  In  addition,  the  strict  constraint  imposed  by  the  resonance 
criterion  permits  the  localization  of  the  excitation  region.  Although  the 
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j;cnc'r;il  idea  of  beat  excitation  can  be  applied  to  various  plasma  inodes,*"^ 
the  majority  of  previous  theoretical  studies  have  considered  the  generation  of 
electrostatic  modes  in  plasmas  which  are  transparent  to  the  pump  waves.  Here 
we  consider  a  problem  in  which  the  existence  of  the  plasma  nonuniformity  is 
advantageously  used  to  localize  the  generation  of  an  electromagnetic  wave  in 
the  whistler  regime. 

The  geometry  of  the  problem  is  sketched  in  Figure  1.  Two  high  frequency 

electromagnetic  waves,  at  closely  spaced  frequencies  and  are  launched 

from  vacuum  into  a  magnetized  plasma  along  the  direction  of  the  zero  order 

density  gradient  In  this  work  we  concentrate  primarily  on  the  case  where 

the  static  magnetic  field  ^  points  against  i.e. ,  ||  VN^.  This  choice 

allows  ready  comparison  with  ongoing  experiments  performed  in  the  northern 
(12) 

auroral  ionosphere. ^  The  vacuum  polarization  of  the  electromagnetic  waves 
is  chosen  to  be  left-handed  in  order  to  avoid  excitation  of  unwanted  electro¬ 
static  nonlinear  phenomena  near  the  plasma  resonance,  i.e.,  at  Wp  =  ,  where 

ujp  is  the  local  electron  plasma  frequency  and  j  =  1,2.  Each  electromagnetic 
wave  penetrates  into  the  plasma  up  to  its  left-hand  cut-off  layer,  i.e.,  where 
u)p  =  -  n) ,  with  fi  the  electron  gyrofrequency  and  it  is  assumed  that 

u).  >  n.  At  this  location  each  wave  is  reflected  and  exhibits  an  Airy  function 
J 

standing  wave  pattern.  The  two  Airy  patterns  are  shifted  with  respect  to  each 
other  because  they  are  centered  about  different  reflection  layers.  It  is 
this  shift  of  the  two  patterns  that  produces  a  spatially  localized  inter¬ 
action  pattern  oscillating  at  the  beat  frequency  w  =  with  a  charac¬ 
teristic  k^  spectrum,  where  k^  is  the  wavenumber  along  The  essential 

idea  is  this  excitation  scheme  consists  of  adjusting  the  various  free  para¬ 
meters  such  that  the  k  corresponding  to  a  whistler  at  w.  -  w,  falls  as 
closely  as  possible  to  the  spectral  peak  associated  with  the  beat  of  the 
individual  Airy  patterns.  Physically,  the  idea  is  to  generate  a  matched 
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antenna  inside  the  plasma  near  the  location  defined  by  u  “  =  w.  (u' .  -  fi) . 

P  J 

The  possibility  of  matching  to  a  whistler  is  greatly  facilitated  because 

this  mode  exhibits  a  dispersion  relation  which  is  essentially  constant  over 

the  region  where  the  beat  antenna  is  formed. 

A  semi -quantitative  criterion  for  the  feasibility  of  generating  a  matched 

antenna  in  a  plasma  can  be  obtained  as  follows.  Near  the  reflection  layer  of 

the  electromagnetic  wave  the  Airy  pattern  exhibits  an  effective  wavenumber 
2/3 

^  (.kjL)  /L  where  =  u^/c,  with  c  the  speed  of  light  and  L  the  density 
gradient  scale  length.  Since  w/w^  <<  1,  the  beat  interaction  pattern  has  a 
significant  spectral  amplitude  at  a  fraction  of  k^.  Equating  the  wavenumber 
of  the  whistler  wave  with  and  using  the  whistler  dispersion  relation 

7  1/2 

k^  =  (w/c)  results  in  the  condition 


(k.L)-2/3 


(1) 


2  2 

where  it  is  assumed  that  Uj  "v  Wp  Equation  (1)  indicates  that  in  plasmas 

in  which  the  density  gradient  is  sufficiently  large  it  is  possible,  in  prin¬ 
ciple,  to  use  beat  excitation  to  resonantly  excite  a  low  frequency  whistler, 
i.e,,  w/fi  •<  1.  For  instance,  in  the  ionosphere  k^L  'v  10^,  S1/2tt  'v  1.4  MHz, 
hence  Eq,  (1)  shows  the  beat  coupling  to  waves  in  the  10  kllz  range  is  possible. 


2 

In  a  laboratory  plasma  with  an  attainable  k^L  'v  lO'^  and  |B^|  'v  400G  coup¬ 
ling  to  whistlers  at  u)/2ti  v  so  MHz  is  expected. 

Although  the  most  efficient  coupling  to  a  whistler  mode  is  accomplished 
by  an  oscillating  current  perpendicular  to  i.e.,  magnetic  multipole  ex¬ 
citation,  such  a  current  cannot  be  directly  generated  by  the  beat  of  two  high 
frequency  waves  propagation  along  Therefore,  in  this  geometry  coupling 

V 

is  accomplished  through  the  less  efficient  electric  dipole  excitation  due  to 
the  generation  of  a  beat  electric  field  along  B  .  The  aim  of  this  study  is  to 
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calculate  the  coupling  coefficient  for  this  process  taking  into  consideration 
the  important  effects  of  plasma  nonuniformity. 

The  beat  electric  field  along  ^  arises  due  to  the  charge  separation 
produced  by  the  beat  ponderomotive  force  acting  on  the  plasma  electrons. 

The  weak  but  finite  variation  of  the  high  frequency  pumps  across  i.e., 
finite  spot  size  at  the  reflection  layer,  provides  the  necessary  coupling 
to  a  slightly  oblique  whistler.  The  transverse  wavenumber  spectrum  of  the 
beam  spot  at  the  reflection  layer  defines  the  spectrum  of  oblique  whistlers 
radiated  by  the  beat  antenna. 

The  theoretical  description  of  whistler  mode  excitation  by  the  beat 
of  two  high  frequency  pumps  requires  three  separate  calculations:  1)  a  suit¬ 
able  representation  of  whistler  mode  properties,  2)  calculation  of  nonlinear 
beat  coupling,  i.e.,  the  source  term  for  the  whistler  wave  equation,  5) 
radiation  of  whistler  using  appropriate  source  functions  in  a  nonuniform 
plasma.  Next  we  proceed  to  discuss  each  of  these  calculations  in  Section 
11.  A  discussion  of  the  results  is  given  in  Section  111. 
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l  I .  I  hoory 

I'irst  we  proceed  to  describe  the  whistler  mode  generated  by  the  beat 

excitation  process.  A  simplifying  l)ut  realistic  assumption  is  that  the 

dispersion  properties  of  the  whistler  mode  can  be  obtained  by  assuming  that 

the  plasma  density  is  constant  and  equal  to  the  density  n^  at  the  reflection 

la>’er  of  the  high  frequency  waves.  This  is  justifiable  because  the  whistler 

mode  is  far  from  any  cut-off  or  resonance,  and  for  the  plasmas  of  interest 

k  1.  ■  1  wliile  the  beat  frequency  antenna  is  effective  over  spatial  scales 

of  the  order  of  The  whistler  wave  at  frequency  u  propagates  at  a 

small  angle  relative  to  B  with  a  wavenumber  k  =  k  +  k  z,  where  k  = 

-o  -  -j_  Z-'  -j_ 

k^  i  X,  fixed  by  the 

finite  transverse  dimension  of  the  high  frequency  beams,  k^  is  determined 
self-consistent ly  from  the  plasma  dispersion  relation  at  frequency  w. 

Using  the  cold  plasma  fluid  equations  together  with  Maxwell's  equations 
yields  the  following  differential  equations  which  determine  the  spatial  depen¬ 
dence  of  the  whistler  electric  field 


-  T  =  k 


k  ±ik 


-r  -  [Ck  ±ikj/2kn  .  [2) 


(k,/k  3')E  =  [(.k  +ik  yk‘l  ^  -  ((k  -ik  )/k^]  ^  , 
o  ^  z  '  V  X  o'  dz  ‘  y  X  o'  dz  ’ 


where  li’(z)  =  (!i  +  il;  )/2,  l  =  1  -  (m  /u))-  [1  ±(il/u'))  ^  ,  c,,  =  1  -  (u)  /u))", 

X  \  I  ^ 

k  =  ./c,  I:  and  b  arc  the  x  and  y  components  of  the  electric  field,  w  is 
o  X  y  po 

the  electron  plasma  frequency  at  the  reflection  layer,  and  a  time  dependence 
cxp(-iait)  is  assumed.  The  plus  and  minus  superscripts  refer  to  the  right-  and 
left-handed  polarizations,  respectively.  The  z  dependence  is  retained  in  F.qs. 
(2)  and  (3J  in  order  to  incorporate  the  dipole  field  generated  by  the  beat 
frequency  antenna.  For  small  angles  of  propagation  relative  to  the  magnetic 


-7- 


ficld  lines,  one  obtains  from  Eq.  (2)  the  approximate  whistler  dispersion 
relation 


(k  /k  ^  ^  ^  w 

z  o  po 

which  is  seen  to  be  insensitive  to  kj^.  The  quantity  is  much  larger 

than  unity,  and  the  propagating  mode  in  Eq.  (4)  is  associated  with  the  F." 
wave  while  the  evanescent  mode  corresponds  to  the  E^  wave.  The  high  evanes¬ 
cence  of  the  latter  and  its  weak  coupling  to  the  propagating  mode  allow  us  to 
neglect  its  effect  on  the  propagation  of  the  E~  wave.  Incorporating  these 
assumptions,  one  finds  that  the  whistler  mode  behavior  is  governed  by  the 
following  scaled  differential  equation. 


.  2_-  „ 

— -  +  b  E  =  S, 
dc2 

■*  -  1  2  2 
where  b"  =  c  -  '^(J'x/k^')  '  uOp^  /ull,  z,  =  k^z, 

dE 

S  =  -[(k  -  ik  )/2k  J  . 

y  X  o^  dc 


(5) 


(6) 


The  source  term  S  arises  from  the  generation  of  an  ambipolar  electric  field 
^  due  to  the  beat  ponderomotive  force  produced  by  the  high  frequency  pumps. 
This  effect  is  calculated  next. 

As  already  stated  in  Sec.  I,  the  relevant  nonlinearity  in  the  problem 
is  the  beat  of  the  two  high  frequency  pumps  near  their  reflection  layers. 

The  lowest  order  nonlinear  effects  occur  at  zero  frequency  (d.c.),  at  low 
frequency  u  =  Wj  -  and  at  high  frequencies  ojj  +0)2*  ^“2’ 

frequency  effect  is  produced  by  the  self-modulation  of  each  wave,  giving 
rise  to  a  modification  of  the  zero  order  density  profile  which  we  do  not 
include  in  the  present  calculation.  The  high  frequency  sideband  at 
as  well  as  the  harmonics  2a)j,  2ui^  are  unable  to  excite  a  wave  because  the 
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spatial  interaction  pattern  does  not  provide  a  suitable  match,  as  is  pos- 
sil)le  for  the  whistler  wave  at  -  (1)2. 

The  nonlinearity  at  lo  =  Wj  -  m,  can  be  obtained  from  the  electron  fluid 
equation  of  motion  by  iteration.  First  one  neglects  the  nonlinear  effects 
and  calculates  the  high  frequency  response  of  the  electrons  to  each  pump 
separately 


V  =  exp(-i(jj^t)  +  v^^  exp(-iio^t),  (7) 

where  v.  =  u-  ‘  E.,  i  =  1,  2,  and 
_j  -y  ' 


with  Yj  =  il/wj ,  and  e  and  m  are  the  electron  charge  and  mass.  The  nonlinear 
effect  is  calculated  by  inserting  the  linearized  velocities  into  the  averaged 
force  equation.  This  yields  a  ponderomotive  force  oscillating  at  frequency  u> 
acting  on  the  electrons  and  given  by 


F  =  -  n^[m(Vj 


-r-  v-  +  V 
z  dz  — 2 


2z  "  I2  " 


where  B.  is  the  magnetic  field  of  wave  j.  We  choose  to  excite  both  pump 

waves  with  left-handed  circular  polarization  to  avoid  additional  effects 

at  =  Wj .  Tlie  appropriate  representation  of  the  pump  electric  field  is 

then  Ej ,  defined  in  a  manner  similar  to  the  field  of  the  whistler  wave. 

Using  Faraday's  law  to  eliminate  yields 
2 


F  =  - 


2n  e  ,  ,  *  , 

o  NT  1  r.-  d  .  1 

nuo 


jU.2^^(1-Yj)  ^1  dz  ^2  ^  TKYp' *^2  dz  Ei]exp(-ia)t) 


(10) 


where  it  is  seen  that  the  force  points  along  the  density  gradient  VN  and 

o 


that  both  jiumps  have  symmetrical  contribution,  as  expected.  The  effect  of 
this  force  is  to  produce  an  ambipolar  electric  field  oscillating  at  w  whose 
self-consistent  strength  is  found  as  follows. 

The  low  frequency  electron  equation  of  motion  takes  the  form 


m  =  -e  E  +  F/n  , 
dt  — s  —  o 


where  is  the  low  frequency  component  of  the  electron  velocity,  and  ^ 
is  the  self-consistent  ambipolar  field  which  also  acts  on  the  ions  of  mass 
M,  charge  q,  and  velocity  through 
du. 

M  =  q  C  ,  (12) 


where  the  effect  of  the  ion  ponderomotive  force  is  neglected  because  of  the 
large  inertia  of  the  ions  and  it  is  also  assumed  that  w  >>  (i.e.,  neg¬ 

ligible  resonance  enhancement  from  the  ion  gyrofrequency  i'2^') .  Solving  for 
the  linearized  velocities  from  Eqs.  (11)  and  (12)  and  using  the  continuity 
equations  together  with  Poisson's  equation  yields 

-(w  /lo)^ 

E  =__EO -  p 

VS  - 

2  2 

where  =  1  -  (lo^^/u)  -  (Up^/w)  ,  is  the  ion  plasma  frequency,  and 

•7 

characteristically  (iOpj^/ui)‘  <<  1. 

It  is  now  possible  to  use  Eqs.  (10)  and  (13)  to  obtain  the  source 


term  S  in  Eq.  (6) 


S(k^,  ky,  c)  = 


moJ.ai-G 

12s 


dT  dr^2  " 


1  / 

I-Y2  '^’2  d^  ^ 


where  the  r  .  ituue  of  each  pump  field  is  expressed  as 
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(15) 


riie  quantity  ^  corresponds  to  tlic  electric  field  amplitude  of  wave  j  at 
at  it^  ret  lection  layer,  and  tlic  function  f(k^,ky)  in  liq.  (14)  is  the  Fourier 
transform  of  the  product  of  the  transverse  envelopes  of  each  pump  beam,  i.e., 
tj(x.y)  in  Fq .  (15).  A  tvqucal  functional  dependence  is  f ^ (x,y)  =  exp[-(x  +y“)/d“] 
with  k  d  ■  1,  and  kj_ 2/d. 

I'o  completely  determine  the  spatial  dependence  of  the  source  term  S 
we  need  to  find  the  functions  'i'j  («’'.),  i.e.,  the  wave  equations  for  the  high 
frequency  pumps  need  to  be  solved  near  their  respective  reflection  layer. 

This  is  easily  done  in  the  relevant  limit  k_L/kj  <<  1,  which  gives 

-2  ,2  (0  2 

k  ^  ^  (z)  +  [1  -  ^  (z)  =  0,  (16) 

•’  dz“  J  u).^(l-Y.)  J 

J  J 

where  the  explicit  spatial  dependence  of  the  density  profile  is  indicated, 
h'e  take  this  variation  to  be  essentially  linear  near  the  reflection  layer, 
i.e.,  we  choose 


1 


..“(1-Y.) 


£ 

L  ’ 


(17) 


with  z  =  0  corresponding  to  the  reflection  layer  of  the  higher  frequency 
pump  Substituting  Eq.  (17)  in  Eq.  (16)  gives  the  Airy  differential 

equation.  After  scaling  the  spatial  variable  in  the  equation  of  pump  2 
to  that  in  the  equation  of  pump  1,  one  obtain^  the  solutions 


(il 

=  Ai(/), 

(18) 

=  Ai(f,  +  A), 

(19) 

where  C  =  (kjL)^^^(2/L) ,  A  =  (kjL)  (l-Y^)  -  (^^/‘^^“(l  -  Y^)],  and  Ai 
refers  to  the  Airy  functions.  These  solutions  are  valid  in  the  neighborhood 
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ot  tlie  cut-off  for  the  left  wave.  Using  Cqs.  (18)  and  (19)  in  fq.  (14) 
determines  the  source  term  completely. 

The  next  task  consists  of  solving  Eq.  (5)  with  the  given  source  S  to 
determine  the  amplitude  of  the  whistler  waves  excited  by  the  beat  frequency 
antenna.  Using  the  Green's  function  method  yields 


(k  -ik  )  fC 

E  '^y’  - - 4k — ~  fexp(ibc)  d^ '  exp(-ib?') 


-exp(-ibc)  d^'  E  (c')exp(ibc')}  , 

J  ^ 


from  which  the  asymptotic  fields  are  found  to  be 


E'(k  ,  k  ,  ^  a(±<«)  exp(±ibc).  c  ->  ±  « 

X  y 


where 


0(1®)  =  ±a  d^  [-j-^AiCC)  Ai’(C+A)  +  Ai  (C+A)Ai '  (C)  ]exp(±iKC)  , 


-(w  /w)  (k  -ik  )f(k  ,k  1  eE„,E„_ 

^  po  ‘  ^  y  x’  y  ,  01  02., 

^  e  2  ^  mu)  u  ’ 

s  12 


K  =  bk^L/(kjL)‘ 


and  the  prime  means  a  derivative  with  respect  to  C.  Recalling  the  assumption 
that  fi/wj  <<  1,  this  implies  that  (1  -  Yj)  2  (1  -  Y2) ,  hence  the  two  terms 
in  the  integrand  of  Eq.  (22)  have  nearly  identical  coefficients  and  this 
allows  us  to  approximate  Eq.  (22)  by 

a(“)  =  (24) 

1 

where 


1=1  d?  Ai(C)  Ai(C+A)exp(-iKC)  . 
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To  obtain  I  is  replaced  by  its  complex  conjugate. 

Since  the  Airy  functions  are  not  square  integrable  in  the  interval 

0).  the  l-ourier  transform  integral  I  is  not  well  defined.  The  diffi¬ 
culty  arises  as  ^  -eo  where  the  function  Ai  oscillates  rapidly  but  its 
amplitude  decays  slowly.  This  beha'dor  results  from  our  usage  of  the  Airy 
function  solutions  away  from  the  neighborhood  of  the  turning  points.  Since 
the  antenna  oscillating  at  the  beat  frequency  is  produced  by  the  interaction 
of  the  main  lobes  of  the  Airy  patterns  within  a  distance  of  a  few  whistler 
wavelengths  (2  the  difficulty  introduced  by  the  highly  oscillatory 

tail  is  easily  overcome  by  introducing  in  the  integrand  of  Eq.  (25)  a  weight 

•>  2 

function  w(4)  =  exp(-C“/4A  )  which  is  nearly  constant  over  the  important  inter¬ 
action  region  and  cancels  the  slowly  decaying  unphysical  behavior  of  the  tail 
as  C  -*■  Other  choices  for  w(C)  have  been  tried  and  the  results  are  found 

to  be  relatively  insensitive  to  the  particular  form  used.  To  illustrate  the 
spatial  dependence  of  the  beat  frequency  antenna  field  generated  by  the 
pump  waves,  a  plot  of  the  function  Ai(C)Ai(C  ♦  A)w(C)  is  shown  in  Fig.  2 
for  a  choice  of  =  10.  The  arrows  at  Wj  and  indicate  the  turning  points 
for  each  pump,  and  the  pure  sinusoidal  wave  at  the  top  of  the  figure  repre¬ 
sents  the  wavelength  of  the  excited  whistler.  The  normalized  spectrum 
of  the  antenna  is  shown  in  Fig.  3,  where  the  arrows  indicate  the  corresponding 
wavenumbers  of  the  excited  whistlers;  only  those  wavenumbers  satisfying  the 
whistler  wave  dispersion  relation  of  Eq.  (4)  are  allowed  to  propagate.  The 
ordinate  represents  the  magnitude  of  the  integral  1  and  the  two  sidebands  cor¬ 
respond  to  the  two  possible  excitation  directions  of  the  whistlers  for  each 
given  value  of  kj^.  Since  it  is  expected  that  the  transverse  beam  spot  is  axially 
symmetric  in  kj_,  the  antenna  radiates  the  whistler  waves  along  a  cone  whose  axis 

is  along  the  z-axis  and  whose  half-angle  is  given  by  tan'^(k,/k  )  as  sketched 

z 
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in  l'ii;uro  1. 

Onco  tho  froi|iioiu'y  of  ono  of  the  pumps  is  chosen,  and  the  radiated  whistler 
•  is  selected  then,  for  a  t’iven  de/isity  scale  lengtli  L,  both  K  and  A  are  uniipicl>' 
doterinined.  This  is  shown  in  l''ij>.  -1  where  both  K  and  A  are  plotted  versus  the 
wliistler  frequency  for  fixed  values  of  u)j/2n  and  L  corresponding  here  to  the 
case  ot  ionospheric  plasmas.  The  process  of  fixing  A  also  determines,  indepen¬ 
dently  ot  K,  the  interaction  pattern  of  Fig.  2  and  its  spectrum  in  Fig.  5.  The 
most  efficient  coupling  occurs  when  the  value  of  K  shown  in  Fig.  4  coincides 
with  tlie  maxima  of  the  spectrum  curve  shown  in  Fig.  5.  Although  there  is  no 
a  priori  reason  for  the  value  of  K  to  be  close  to  the  spectral  peak  of  the 
antenna,  it  turns  out  that  for  parameters  typical  of  ionospheric  and  labora¬ 
tory  plasmas  it  is  possible  to  approach  the  optimum  value,  as  exemplified  by 
the  arrows  in  Figure  3. 

Finally,  we  proceed  to  calculate  the  strength  of  the  radiated  signal  and 
the  efficiency  of  the  coupling  between  the  high  frequency  pump  waves  and  the 
whistler  wave. 

We  substitute  Eq.  (24)  into  Eq.  (21)  and  take  the  inverse  Fourier  trans¬ 
form  of  tlie  transverse  dependence  to  obtain 


ir(x,y,z)  'V  E|^  exp[i(+bk^z  +  41  +  0  +  tt/2)], 


2  -♦■  +  CO 


(26) 


where 


-(.  /.)  20Ki;„,u  HI 
hi  ■  — - H  ‘ 

s  1  2  1  d 

-1, 


(27) 


4  is  the  phase  of  I,  and  0  =  tan"  (y/x) .  Equation  (27)  gives  the  total 
amplitude  of  the  radiated  field  ol'  the  whistler  wave.  A  qualitative  sketch 
is  sliown  in  Fig.  5.  It  exhibits  a  null  on  the  2  axis,  and  a  maximum  along 
a  circular  ring  of  radius  d/2. 
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riu'  average  energy  flux  radiated  in  either 

1  -  -* 

lirect  u>n  is  g.iven  in  MkS  units  hy  B 

•\[>i'ess i  on  for  the  total  power  radiated 


the  ])ositive  or  negative  z 
),  whicli  leads  to  the  following 


=  2l)f^  / p 
o  o 


dxdy  n 


R 


2nbi'c  /p 
o  o 


^‘■'^*'01  ^'02  I 


I  /2mojjU)2(l 


-Yj)],' 


where  and  p^  arc  the  permittivity  and  permeability  of  free  space,  respect¬ 
ively  . 


Assuming  that  the  pump  waves  originate  from  a  transmitter  in  vacuum, 


one  can  relate  the  amplitude  to  the  vacuum  pump  power  through  the 

2  _  2  2 

approximate  expression  ^  Tid  ^^01  *  ^02  ^  efficiency  with 

which  the  whistler  wave  is  excited  by  the  pumps  is 


n  =  P^/Pp.  C29) 

This  is  a  nonlinear  efficiency  because  it  depends  on  the  pump  power.  In  order  to 
obtain  a  more  practical  expression  for  n  we  rewrite  Eq .  (29)  as 

(u,  /4;n)(eKli|)^ 

n  =  - 4-rbl - T"  P  .  C30) 

Slid  i^/Pq  [mw^w^ci-vp]"  P 


where  we  made  the  reasonable  assumption  that  E^^  =  E^^  =  Equation  (30) 

gives  the  coupling  efficiency  ti  as  a  function  of  the  total  power  radiated 
by  the  pump. 


1 
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1  1 1  .  I)  i  SL  uss  i  on  \ 

Ihc  present  calculation  provides  a  clear  example  of  how  the  nonlinear 
efficiency  for  the  radiation  of  a  low  frequency  electromagnetic  signal  due 
to  the  beat  of  two  high  frequency  pumps  can  be  determined  from  first  princi¬ 
ples.  The  calculation  includes  plasma  nonuniformity  as  well  as  finite  spatial 
extent  of  the  pumps.  The  figure  of  merit  for  the  beat  process  is  given  by  the 
compact  expression,  shown  in  Eq.  f30).  A  rough  order  of  magnitude  evaluation 
of  this  expression  shows  a  radiated  power  whose  dependence  on  pump  power 

P  scales  as  "v  (v  /c)^P  ,  where  v  refers  to  the  oscillating  velocity  of 
p  R  ^  os  p  os  b  j 

an  electron  due  to  the  high  frequency  pump,  and  c  is  the  speed  of  light. 

This  dependence  illustrates  readily  the  extreme  inefficiency  of  the  process 
at  low  pump  power  (i.e.,  small  .  Another  fundamental  factor  which  lowers 
the  efficiency  is  that  the  coupling  to  the  whistler  wave  occurs  through  elec¬ 
tric  dipole  excitation.  The  dependence  due  to  this  effect  appears  through 
-4 

the  d  term  in  Eq.  (30). 

A  more  efficient  excitation  mechanism  would  be  to  use  magnetic  coupling, 
i.e.,  to  generate  a  nonlinear  transverse  current  oscillating  at  the  beat  fre¬ 
quency  u).  However,  within  the  constraints  of  the  present  geometry  (i.e., 

I  I  VN^)  this  more  desireable  situation  cannot  be  realized  directly. 

The  most  advantageous  feature  of  the  beat  excitation  described  here 
is  the  controllability  of  the  spatial  location  from  which  the  whistler  eman¬ 
ates.  By  proper  adjustment  of  the  pump  frequency  Wj  and  usage  of  Fig.  4  in 
conjunction  with  Fig.  3  the  optimum  efficiency  can  be  realized  for  a  desired 
remote  location  in  the  plasma.  Due  to  the  dispersive  properties  of  whistlers, 
the  radiated  low  frequency  signal  exhibits  a  characteristic  headlight  pattern, 
as  is  illustrated  in  Fig,  5,  The  signal  is  generated  at  the  reflection  layer 


of  the  high  frequency  pumps  and  is  radiated  along  a  hollow  cylinder  concentric 
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with  the  pump  beam  illumination  spot. 

Considering  parameters  typical  of  the  auroral  ionosphere,  L  =  50  km,  n/2ii  = 
1.4  Mil;,  reflection  layer  at  150  km,  and  present  transmitter  capabilities,  i.e., 
1’^^  =  5  MW,  =  5  Mil;,  u)^/2ti  ~  4.9  MHz,  and  antenna  gain  of  30  db,  results 

in  values  of  the  scaled  quantities  K  =  4.64,  A  =  10.3,  and  |l|  =  1.4.  These 
conditions  give  rise  to  the  radiation  of  a  whistler  at  ci)/2n  =  100  KHz  with  an 
electric  field  amplitude  of  5yV/m  at  an  extremely  low  nonlinear  efficiency  level 
of  n  ''  lO"^^.  In  spite  of  the  low  system  efficiency,  present  experiments  detect 
signals  at  comparable  levels. 

It  should  also  be  mentioned  that  naturally  occuring  atmospheric  lightning 
radiates  a  broad  spectrum  of  electromagnetic  waves.  In  principle,  the  high 
frequency  noise  waves  can  trigger  the  secondary  radiation  of  whistlers  by  the 
mechanism  described  in  this  work. 

For  parameters  realizable  in  a  large  laboratory  plasma  device,  |l^|  = 

400C.,  u)j/2n  =  2.8  GHz,  u)^/2ii  =  2.5GHz,  L  =  10m,  one  obtains  K  =  4.33,  A  =  11.2, 

|li  -  1,  which  results  in  the  excitation  of  a  whistler  at  w/2n  =  300  MHz  with 

-8 

an  electric  field  amplitude  of  50  mV/m  at  a  coupling  efficiency  of  n  10  for 
an  input  pump  power  of  1  KW.  Thus,  the  effect  is  amenable  to  experimental 
investigation. 

Finally,  it  should  be  noted  that  if  the  density  gradient  points  at  a 

small  angle  0  relative  to  the  magnetic  field,  then  the  results  described  in 

2  2 

this  stud)  remain  valid  up  to  an  order  (fi/uij)  sin  6.  For  ionospheric  plasmas 
this  condition  restricts  the  applicability  of  the  present  analysis  to  high 
latitude  auroral  plasmas. 
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lUgurc  Captions 


I  ij;.  I  iiconictry  of  the  problem,  where  A  is  the  scaled  separation  of  the 
pumps  obtained  in  Eq.  (19). 

Fig.  2  Spatial  pattern  of  the  beat  frequency  antenna  generated  by  the  two 
high  frequency  pumps.  The  arrows  show  the  turning  points  for  each 
pump  which  are  separated  by  a  scaled  distance  A  =  10.  The  radiated 
whistler  wavelength  is  shown  on  top  for  reference. 

Fig.  3  Normalized  spectrum  of  the  antenna  pattern  shown  in  Fig.  2. 

The  arrows  indicate  the  magnitude  of  I  corresponding  to  the  wave- 
numbers  of  the  excited  whistlers. 

Fig.  4  Normalized  whistler  wavenumber  K  and  separation  of  the  pumps  A 

versus  whistler  frequency  w  for  pump  frequency  u^/2v  =  5  MHZ,  and 
density  scalelength  L  =  50  km. 

Fig.  5  Sketch  of  the  radiated  whistlers  showing  the  transverse  dependence 
of  the  amplitude  [p  =  (x^  + 
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